incomplete_gamma_upper Function

public pure elemental function incomplete_gamma_upper(a, x) result(rst)

Computes the upper incomplete gamma function.

The upper incomplete gamma function is defined as: $$\Gamma(a, x) = \int_{x}^{\infty} t^{a-1} e^{-t} \,dt$$

Remarks

The value overflows for arguments beyond roughly $a = 171$, where $\Gamma(a)$ itself exceeds the range of a double precision number. Use regularized_gamma_upper where only the ratio to $\Gamma(a)$ is required.

See Also

• Wikipedia

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=real64),	intent(in)			::	a	The coefficient value.
real(kind=real64),	intent(in)			::	x	The value at which to evaluate the function.

Return Value real(kind=real64)

The function value.